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Abstract. While the implementation of single particle coherent diffraction imaging 
for non-crystalline particles is complicated by current limitations in photon flux, hit 
rate, and sample delivery a concept of many-particle coherent diffraction imaging offers 
an alternative way to overcome these difficulties. Here we present a direct, non-iterative 
approach for the recovery of the diffraction pattern corresponding to a single particle 
using coherent x-ray data collected from a two-dimensional (2D) disordered system of 
identical particles, that does not require a priori information about the particles and 
can be applied to a general case of particles without symmetry. The reconstructed 
single particle diffraction pattern can be directly used in common iterative phase 
retrieval algorithms to recover the structure of the particle. 
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1. Introduction 

It was recently realized [H El E] that analysis of diffraction patterns from disordered 
systems based on intensity cross-correlation functions (CCFs) can provide information 
on the local symmetry of these systems. A variety of disordered systems, for example, 
colloids and molecules in solution, liquids, or atomic clusters in the gas phase can be 
studied by this approach. It becomes especially attractive with the availability of x- 
ray free-electron lasers (FELs) [U El E]. Essentially, scattering from an ensemble of 
identical particles can provide, in principle, the same structural information as single 
particle coherent imaging experiments [7J that are presently limited in the achievable 
resolution [8] . Clearly, a large number of particles in the native environment can scatter 
up to a higher resolution, at the same photon fluences, comparing to experiments on 
single particles injected into the FEL beam. However, it is still a challenge to recover 
the structure of individual particles composing a system using the CCF formalism. 
In the pioneering work of Kam [HI HO] it was proposed to determine the structure of 
a single particle using scattered intensity from many identical particles in solution. 
However, this approach, based on spherical harmonics expansion of the scattered 
amplitudes, was not fully explored until now. Recently, it has been revised theoretically 
[XT] IT2"1 [13] and experimentally [14] . The possibility to recover the structure of individual 
particle was demonstrated in systems in two-dimensions (2D) [TTl [121 [14] and in three- 
dimensions (3D) [13] using additional a priori knowledge on the symmetry of the 
particles. Unfortunately, these approaches, based on optimization routines [HI H3] and 
iterative techniques [121 [14] , do not guarantee the uniqueness of the recovered structure 
unless strong constraints are applied. This all leads to a high demand in finding direct, 
non-iterative approaches for recovering the structure of an individual particle in the 
system. 

In this paper we further develop Kam's ideas and propose an approach enabling 
unambiguous reconstruction of single-particle intensity distributions using algebraic 
formalism of measured two- and three-point CCFs without additional constraints. Once 
the single-particle intensity is obtained, conventional phase retrieval algorithms [151 HB] 
can recover the projected electron density of a single particle. Our approach is developed 
for 2D systems of particles which is of particular interest for studies of membrane 
proteins. It can be also used to study 3D systems provided that particles can be aligned 
with respect to the incoming x-ray beam direction. 

2. Basic equations 

We consider a scattering experiment in transmission geometry [see figure [I] , where the 
direction of the incoming x-ray beam is perpendicular to the 2D sample plane, and 
simulate a set of diffraction patterns. Our sample consists of an arbitrary small number 
N of identical, spatially disordered particles [see figure |2]( a)]. Also, we assume a uniform 
distribution of orientations of particles in the system. First, for simplicity, we consider 
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a model where the total scattered intensity I(q) is represented as an incoherent sum of 
intensities 1^. (q) corresponding to individual particles in the system, 

N 

= (i) 

i=l 

where q is the momentum transfer vector and ipi is the orientation of the i-th particle. 
This model is a good approximation for dilute systems when the mean distance between 
particles is much larger than their size [21 E] • In this approximation we neglect the inter- 
particle correlations due to coherent interference of scattered amplitudes from individual 
particles. Later, in simulations, we generalize our approach to the case of coherent 
scattering from a system of particles in the presence of Poisson noise and demonstrate 
the applicability of our approach. 



Incident 
x-ray beam 
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Figure 1. (Color online) Geometry of the diffraction experiment. The incident x-ray 
beam coherently illuminates a 2D disordered sample and produces a diffraction pattern 
on a detector. The direction of the incident beam is defined along the z axis of the 
coordinate system. 

In the frame of kinematical scattering, the intensity (q) scattered from a single 
particle in some reference orientation ipo is related to the electron density of the particle 
PVo( r ) by the following relation 



p^ (r)exp(iqrr)dr 



(2) 



Once the single-particle intensity I^ {q) is determined, conventional phase retrieval 
algorithms [151 EE] can recover the projected electron density p^ (r) of a single particle. 
Our goal is to determine the scattering pattern of a single particle I^ {q) using a large 
number of diffraction patterns I(q) corresponding to different realizations of the system. 

The intensity (q) scattered from a single particle can also be described in a polar 
coordinate system q = (q, if) as a function of the angular position (0 < if < 2tt) on 
the resolution ring q [see figure p^b)]. For each resolution ring q, one can represent the 
intensity function (q, if) as a Fourier series expansion, 



(3) 
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Figure 2. (Color online) (a) 2D disordered system composed of N = 10 clusters in 
random positions and orientations (b) Instantaneous scattered intensity I(q, ip) as a 
function of the angular position <f> on the resolution ring q. (c),(d) Definition of the 
momentum transfer vectors in the derivation of the two-point C(qi, q 2 , A) (c), and the 
three-point CCFs C(q 1 ,q 2 ,q 3 ,A 1 ,A 2 ) (d) . 



where 7™^ o are the Fourier components of I^ (q,(p). 

For a 2D system of identical particles, the intensity Iy (q,ip) scattered from a 
particle in the reference orientation i/; — is related to the intensity 1^. (q, ip) scattered 
from a particle in an arbitrary orientation ipi as I^(q, ip) = Iip (q, ip — ipi)- Applying the 
shift theorem for the Fourier transforms [17] we obtain for the corresponding Fourier 
components of the intensities, J^. = I™^, exp(— inipi). Using these relations we can 
write for the Fourier components I™ of the intensity I(q, ip) scattered from N particles 

N 

in q = Cfc E exp(-in^) = I^ A n , (4) 
i=i 

where A n = J2iLi exp(— inipi) is a random phasor sum [18]. According to ^ the 
intensity scattered from a single particle can be uniquely determined by the set of 
complex coefficients {Iq^ } = {\Iq^ \, 4>q^ Q = arg(/^ o )}. Here we propose a direct 
approach for determination of these Fourier components {Ig^ } applying two- and three- 
point CCFs to the measured intensities I(q, ip) scattered form N particles. 

We start with the two-point CCF defined at two resolution rings q% and q<i [§1 El [2] 

C( qi , q 2 , A) = (J( qi ,<p)T(q 2 , <p + A)) v , (5) 

where < A < 2tt is the angular coordinate [see figure [2^c)], I(q,ip) = I(q,ip) — 
(I(q,ip)) is the intensity fluctuation function, and (. . .) denotes the average over the 
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angle if. It can be directly shown [2] that the Fourier components C" q2 of the CCF 
C(qi, q2, A) for n ^ are defined by the Fourier components I™ of the intensities I(q, ipj^ 

cs* = ^r-c ( 6 ) 

and for n = Fourier components (7? g2 = according to the definition of the intensity 
fluctuation function I(q, if). Using M) in (pi) and introducing statistical averaging (. . .) M 



over a large number M of diffraction patterns one can get 

|2\ 



fin \ _ jn* jn . / I A | 2 \ - T n * T n . AT (7) 

Here, we used the fact that for a uniform distribution of orientations of iV particles 
(|A„| 2 ) M asymptotically converges to N for a sufficiently large number M of diffraction 
pattern^]]] 

Equation (J7| can be used to determine both, the amplitudes \Iq^ Q \ and phases 
(fig^o (f° r n. > 0) of the Fourier components J"^ o associated with a single particle. For 
example, to determine the amplitudes equation (j7]) can be applied successively to three 
different resolution rings q±, q2 and ^3 connecting each time a pair of q- values. Direct 
evaluation gives for the amplitudes of the Fourier components of a single particle on the 
ring qi 

IOol=W^ ( 8 ) 



where X™^ 



91.92/ m 



^93,91) I / ( ^93.92) is an experimentally determined 
quantity. Applying (|8j) to different resolution rings q and orders n all required 
amplitudes |/^ | can be determineoj^J Equation ([8| should be used with care to avoid 
possible instabilities due to division by zero. For this purpose one should exclude from 
consideration the cases when (C^ | is close to zero. Since the Fourier components 
of the intensities obey the symmetry condition 7™^, o = I~^ Q it is sufficient to determine 
Iqip for n > 0. The 0-th order Fourier component by its definition [2j [3] is a real- valued 
quantity, J° = (I(q,if)) , and can be determined from the experiment as well. Using 
Q one can readily find that ij^ = Jj/JV. 

Equation ([7]) also determines the phase difference between two Fourier components 
Iq^g and Ig\^ of the same order n, defined at two different resolution rings q\ and q2, 



= ^o-C^o- (9) 
Notice, that one can freely assign an arbitrary phase to one of the Fourier components 
Iqjfg, which corresponds to an arbitrary initial angular orientation of a particle. 

§ We note that the Fourier components of the intensity fluctuation function 7™ can be expressed 
through the Fourier components of intensity I™ using relation I™ = I™ — I® ■ d n .o, where 5 n $ is the 
Kronecker symbol. 

|| Note, that in this paper, contrary to [2J|3], we use a not normalized CCFs which, in particular, results 
in different asymptotic values of (|A n | 2 ) M . 

If number of particles in the system is not known a priori then N in ^ has to be considered as 
a scaling factor that is obtained on the final stage of reconstruction of a single particle intensity (see 
section |3l . 
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Assuming, for example, <j>^ u ^, = and using (|9j) one can directly determine the phases 
of the Fourier components with the same n-value on all other resolution rings q 2 7^ q\ 

To completely solve the phase problem, it is required to obtain additional phase 
relations between Fourier components with different n values on different resolution 
rings q. These relations can be determined using a three-point CCF introduced by Kam 
P HO]- The important aspect of our approach is to use the three-point CCF defined on 
three different resolution rings [see figure |2](d)], 

C(q u q 2 , q 3 , A 1; A 2 ) = (%, <p)I(q 2 , <p + A 1 )I(q 3 , if + A 2 )\ , (10) 



<p 

contrary to [11] [12] where the three-point CCF was defined on two resolution rings. 
Similar to C™ q2 it can be shown (see Appendix) that for the Fourier components of this 
CCF the following relation is valid, 

rini,ri2 _ r(ni+ra 2 )* rni rn 2 

9l>92,93 91 12 93 ' V / 

for n% 7^ 0,?22 7^ 0,ni 7^ — n 2 . Using ^ in (11) and performing statistical averaging we 
obtain for the averaged Fourier components of the three-point CCF 

/ /~tni,ri2 \ _ Ani+m)* rni rn 2 l\ \ /iq\ 

where A nii „ 2 = X^f^-,fc=i exp{i[(ni + r^)^ — niipj ~ ^ipk)}- Our analysis shows (see 
Appendix) that the statistical average (A ni n2 ) M converges to N for a sufficiently large 
number M of diffraction patterns, i.e. arg[(A nii „ 2 ) M ] = 0, and we get from (12) the 
following phase relation 

«K[(C£C) M ] = Co + Co - <CT 2) - (13) 



Equation (13) determines the phase shift between three Fourier components lj: ni 1 ""' 



91,^0 ' 

I&ifo> anc ^ Iqaifo °^ different order n defined on three resolution rings. If n\ = n 2 = n and 
n 3 = 2n, equation (13) reduces to a particular form, giving the phase relation between 
Fourier components of only two different orders n and 2n, 

a^tea.) J = + K,+> ~ ( 14 ) 
Equations ([8]), ([9]) and (13) constitute the core of our approach and allow us to 
directly and unambiguously determine the complex Fourier components J^ o using 
measured x-ray data from a disordered system of N particles. The obtained Fourier 
components J^ can be used in (|3| to recover the scattered intensity I^, (q,(p) 
corresponding to a single particle. 



3. Recovery of the projected electron density of a single particle 

We demonstrate our approach for the case of a coherent illumination of a disordered 
system of particles, in the presence of Poisson noise in the scattered signal. We recover 
diffraction patterns and projected electron densities for two different particles, a centered 
pentagonal cluster [figure [3t a)] that has 5-fold rotational symmetry and an asymmetric 



(g)l (h) 
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Figure 3. (Color online) (a),(b) Scattered intensity (logarithmic scale) calculated for 
a single pentagonal cluster (a) and an asymmetric cluster (b) (clusters are shown in the 
insets). (c),(d) Coherently scattered intensity from a disordered system consisting of 
N = 10 clusters in random position and orientation. Scattered signal corresponding to 
the incident fluence of 10 12 and 10 13 photons/25 /mi 2 for pentagonal and asymmetric 
clusters correspondingly. Diffraction pattern shown in (d) also contains Poisson noise. 

(e) ,(f) Scattered intensity corresponding to a single pentagonal (e) and asymmetric 

(f) clusters recovered from M = 10 5 diffraction patterns of the form (c) and (d) 
correspondingly. (g),(h) Structure of a single cluster reconstructed by an iterative 
phase retrieval algorithm using the diffraction patterns shown in (e) and (f). 
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cluster [figure p[b)]. Both clusters have a size of d — 300 nm and are composed of 
polymethylmethacrylate (PMMA) spheres of 50 nm radius. 

In both cases we consider kinematical coherent scattering of x-rays with wavelength 
A = 1 A from a system of iV = 10 clusters in random positions and orientations, 
distributed within a sample area of 5 x 5 jjm 2 [see figure [2^a)]. Diffraction patterns are 
simulated for a 2D detector of size D = 24 mm (with pixel size p = 80 yum), positioned 
in the transmission geometry at L = 3 m distance from the sample [see figure [I] . This 
experimental geometry corresponds to scattering to a maximum resolution of 0.25 nm -1 . 
For given experimental conditions the speckle size corresponding to the illuminated area 
is below the pixel size of the detector. At the same time the speckle size corresponding 
to the size of a single particle is about 12 pixels, that provides sufficient sampling for 
phase the retrieval. 

The coherently scattered intensities simulated for single realizations of the system^ 
are shown in figures [3^c) and|3^d). We note, that in simulations the intensity /(<?, y?) 
scattered from a system of iV particles was calculated as a coherent sum of the 
scattered amplitudes A^. (q, tp) from each particle, i.e, I(q,(p) = \ J2iLi L P)\ 2 - The 

incident fluence was considered to be 10 12 and 10 13 photons/25 /im 2 for pentagonal and 
asymmetric clusters, respectively. In the case of asymmetric clusters additional Poisson 
noise was included in the simulations. The Fourier components of two-point and three- 
point CCFs [equations ^ and (12)] were averaged over M = 10 5 diffraction pattern^ 
Once all amplitudes 



"9,^0 



and phases 



1^0 



of the Fourier components I™ ^ are 



determined, the diffraction pattern corresponding to a single particle can be recovered 
by performing the Fourier transform (|3). Using experimentally accessible quantities 



fl M 



M?,*, ^d X^ 



jn 

q,ipo 



we can rewrite (3) as 



f/ M 



N 



+ E 

71 = — OO 

n^0 



q,ipo 



exp(i#^)exp(irap). 



(15) 



Here we discuss determination of the amplitudes and phases of the Fourier 
components 7™^ o associated with a single particle, using x-ray data from 2D system 
of pentagonal clusters [figure]^ a)] described above. Applying expression Q to different 
q we find all required amplitudes y/N\I^ \ (for n > 0) scaled by the factor y/N. The 
0-th order Fourier component is determined as J°^ o = ((I{q, f))^) M /N. In figure jijthe 
values \/N\I™y Q \ derived from (8) normalized by \(I(q, ( f)} ip / are shown for n < 40 at 
three different resolution rings q\ = 0.21 nm -1 , q 2 = 0.23 nm -1 and = 0.25 nm -1 . 
Due to 5-fold symmetry of the pentagonal cluster, in the g-range accessible in our model 

+ All simulations of diffraction patterns were performed using the computer code MOLTRANS. 

* According to our simulations, two-point and three-point CCFs defined on different resolution rings 



q%, q2 and (73 [see |5| and (10)] even in the case of coherent illumination of a disordered system do not 
contain the inter-particle contribution. Contrary to that, as it was shown in [51 [3], this is not the case 
for two-point CCFs defined on the same resolution ring q\ = q2, when inter-particle contributions and 
noise can give a substantial contribution. Therefore, in our present work we avoided this last case for 
the unique determination of amplitudes and phases of the Fourier coefficients of single particles. 
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experiment we need to determine the phases of the Fourier components I"w, Q only with 
n — 10 • I, 1 < I < 4, where I is an integer. Below we present the algorithm for the 
phase determination using equations Q, (13) and (14). 
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Figure 4. (Color online) Normalized amplitudes \Ig^ \norm = l^vJAI-^J^) 
determined using Q at three different resolution rings (a) qi — 0.21 nm -1 , (b) 
q2 = 0.23 nm -1 and (c) — 0.25 nm" 1 . 



(i) First, we assign an arbitrary phase to one of the Fourier components Iq^, which 
corresponds to an arbitrary initial angular orientation of a particle. It is convenient 
to assign this phase to the Fourier component with the smallest available n value. 
In our case we assign a zero phase <fql$ — to the Fourier component I^^, Q - 
ii) The phases of the Fourier components with n = 10 on other resolution rings 
Qj U = 1,2) were obtained by using (|9|), i.e., = + arg^C^^ 



(iii) The phase of the Fourier component I^j, on the resolution ring q± was determined 
by using Q with n = 10, <^ = ^ + io^ _ aig^Cj^) J. 

(iv) The phases of the Fourier components with n = 20 on other resolution rings 
Qj U = 1,2) were obtained from (9) <f>f jM = ^20^ + ^[^720^ ^ 

(v) The phase of the Fourier component li ^,. on the resolution ring q\ was determined 



by using (13) with m = 10 and n 2 = 20, <j>™ = #° + #° _ arg [(^ 10 -i 



91,00 ~~ ^92,l/ , ^93,^0 _ cll SL\ , - v l}l,(J2,93 / jyyJ 

The process was continued in a similar way until all phases at each g-value were 
determined. The same procedure has been applied for the case of 2D system of 
asymmetric clusters [figure [3^b)]. Due to asymmetric structure of particles, in the 
same g-range we need to determine the phases of a significantly larger set of the Fourier 
components with n = 2 ■ k, 1 < k < 24, where k is an integer. Assigning a zero 
phase <\? ^ o = to the Fourier component with n = 2 we successively determine the 
phases of the Fourier components of the higher orders up to n = 48. We note, that 
the data redundancy intrinsic to (^v, g2 ) M and (pq^q£ q ^J offers a lot of flexibility in 
determination the possible ways of solving the phases of the Fourier components Ig^ - 

If the number of particles N in the system is not known [see equation ( 15 )] , it can be 
determined using the positivity constraint (q, y?) > for the recovered intensity. One 
can slightly relax this constraint and allow a small fraction of pixels with negative values 
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in order to get an image with a better contrast. These negative values can be substituted 
by zeros before using this diffraction pattern in the iterative phase retrieval algorithms 
for the reconstruction of the particle structure. This procedure can be justified by 
inaccuracies which arise due to statistical estimate of the CCF and also different errors 
in calculations. For example, for the case of pentagonal clusters in our simulations we 
applied the value of N = 9.8, which corresponds to 0.27% of negative pixes on the 
detector [see figure [5] . This value is very close to the theoretically predicted N = 10 for 
a uniform distribution of orientations of particles. 




Scaling factor, TV 

Figure 5. (Color online) Dependence of the relative number of pixels with negative 
values (normalized to the total number of pixels on the detector) in the recovered 
diffraction pattern I^ (q,ip) as a function of the scaling factor N. The red point 
indicates the value of N = 9.8 used in the final reconstruction of intensity I^ (q, ip). 



4. Results and discussion 

The diffraction patterns corresponding to single particles recovered by our approach are 
presented in figures [3]^e) and[3^f). As one can see from these figures, the recovered single 
particle diffraction patterns reproduce the diffraction patterns of individual clusters 
shown in figures [3^a) and ^b) very well. The structure of the clusters reconstructed 
from these recovered diffraction patterns was obtained by a standard phase retrieval 
approach and is presented in figures ^g) and[3^h). In our reconstructions we used an 
alternative sequence of hybrid input-output (HIO) and error-reduction (ER) algorithms 
[15] . The comparison of the single cluster structures obtained by our approach 
[figures |3tg) and ^h)] with the initial model shown in the insets of figures |3](a) and 
|3^b) confirms the correctness of our reconstruction. Both clusters were reconstructed 
with a resolution up to 25 nm. These results clearly demonstrate the ability of our 
approach to recover the single-particle structure from noisy data obtained in coherent 
x-ray scattering experiments. 

The ability of the presented approach to recover the diffraction pattern of a single 
particle relays on the accuracy of the two- and three- point CCFs determination. The 
statistical properties of the CCFs and their convergence to the average values strongly 
depend on the number of particles N in the system, their density and the distribution of 
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their orientations [2, 3J. In particular, the value of the scaling factor N in the expression 
of the Fourier components of the CCF is defined by the statistical distribution of particles 
in the system. We performed a few simulations with varying number N of particles in 
the system: a) for a Gaussian distribution of the number of particles with an average 
value (N) = 20 and standard deviation a = 2; b) for a Poisson distribution with an 
average number of particles (N) = 30. In both cases the CCFs statistically converge to 
the same average values as for the systems with a fixed number of particles, N = 20 and 
N = 30 (for Gaussian and Poisson distributions, correspondingly). The only difference 
we observed is a slightly slower convergence of CCFs in the case of statistical distribution 
of particles compared to the case with a fixed number N. This means, that in the case 
of systems with a varying number of particles one needs to measure a larger number of 
diffraction patterns. 

Particle density also affects the statistical behaviour of the CCFs. In the limiting 
case of a very dilute system the CCFs calculated for a single diffraction pattern are 
defined by independent structural contributions of N individual particles. In this case 
statistical averaging over diffraction patterns is, in fact, acts on the fluctuating terms 
|A n | 2 in (j7j) and A ni ri2 in (12) and convergence of these terms to their statistical 
estimates determines the number M of diffraction patterns required for averaging [3]. 
In the case of a dense system, the CCFs contain a significant inter-particle contribution 
in the same range of g-values, where the structural contribution of individual particles 
is observed [2j [3]. This additional contribution slows down the convergence of CCFs 
(especially of the three-point CCF) and the number of measured diffraction patterns 
should be increased. In the presence of Poisson noise in the scattered signal one needs 
to accumulate even more diffraction patterns. Additional effects of resolution, solution 
scattering, etc. on the measured CCFs were considered in detail in [19]. Results of 
our simulations with PMMA clusters and asymptotic estimates presented in [3] show 
that, in practice, the number of particles in the system should be less than few dozens. 
This allows to achieve statistical convergence of the CCFs using 10 5 — 10 6 diffraction 
patterns. In the disordered systems considered in this work the particle density is 
characterized by (R)/d ~ 5.0, where (R) is the average inter-particle distance. At these 
conditions it was sufficient to use 10 5 diffraction patterns with Poisson noise to achieve 
convergence of the two and three-point CCFs and perform successful reconstruction of 
the particle structure. In practical applications the convergence of the CCFs can be 
directly controlled as a function of the number M of the diffraction patterns considered 
in the averaging [3]. 



5. Conclusions 



In conclusion, we present here a direct, non-iterative approach for the unambiguous 
recovery of the scattering pattern corresponding to a single particle (molecule, cluster, 
etc) using scattering data from a disordered system of these particles including noise. 
Our simulations demonstrate the successful application of this approach to 2D systems 
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composed of particles with and without rotational symmetry. It can be of particular 
interest for studies of membrane proteins which naturally form 2D systems. Our 
approach can be also applied to 3D systems of particles to recover their projected 
electron density, if a specific alignment of the particles along the direction of the incoming 
x-ray beam can be achieved. We have shown that our approach is robust to noise and 
intensity fluctuations which arise due to coherent interference of waves scattered from 
different particles. We foresee that this method will find a wide application in the studies 
of disordered systems at newly emerging free-electron lasers. 
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Appendix 



The Fourier series expansion of the CCF ( 10 ) can be written as 

oo oo 

C(q 1 ,q 2 ,q 3 ,A 1 ,A 2 ) = E E C^e^e*** m ^ -n 2 , (A.l) 

n-^ — — oo ri2 — — oo 

C 1'Z 3 = (^) a jf r jf*C(g ll ( bl g3,A 1 ,A 2 ) e -^e-^dA 1 dA 2 

= / (n 1+ n 2 )* / n 1/ n 2; ni ^ Q, Tfc ^ 0, n X ^ -n*. (A.2) 

where C^^ qa are the Fourier components of the CCF C(qi, q 2 , (73, A 1; A 2 ). In general, 



(A.2) determines a relation between three different Fourier components of intensity 7" 
of the order m, n 2 and n\ + n 2 , defined on three resolution rings, q\, q 2 and q 3 . 

Considering incoherent scattering from N particles, we can rewrite (A.2) in terms 
of the Fourier components of intensity I™^ Q associated with a single particle, 



/nrni,ri2 
91,92,93/^ 



r (ni+n2)* jm jn 



T ni T n ' z • (A ) 



(A.3) 



where (. . .) M denotes the statistical averaging and we introduce a new random phasor 
sum 



N 

L m,n 2 = E exp{i[(ni + n 2 )^i - mipj - n 2 ^ k }}. 

i,j,k=X 



(A.4) 



To understand the statistical properties of (A ni n2 ) M we split it into three contributions, 

(A ni)n2 ) M = (B nijn2 ) A/ + (C ni ,n 2 ) M C^ n i, n 2) M ' (A-5) 
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where (B r 



'ni,n 2 / > 

subscripts z, j and fc 



(B nijn2 ) 



(D nijn2 ) are defined for different combinations of the 



n 2 /M 



Y exp{i[(m +n 2 )ipi - n x ipj - n 2 ip k ]}) = N, (A.6) 



\C nijn2 ) A y 



=j=k 

E 



M 



■ ■■ + 



E ^ 



+ E ^ 



2 ( ^ cos[ni(^ - ipj)} ) +2(5^ cos[n 2 (^i - ipj)} 



\i>j 



M 



\i>j 



M 



+ 2 / ^ cos[(rii + n 2 ) (ipi - ipj)] 
\i>j i 



M 



(D nijn2 ) 



M 



Y exp{i[(ni + n 2 )*pi 



mi/jj - n 2 ^ k }} 



Now we determine the statistical estimates of the terms (C r 




(A.7) 



(Ai 



'ni,ndM all< ^ (■Dni,n 2 ) A f We 

consider an uniform distribution of orientations of particles ip on the interval (— 7r,7r), 
where two different orientations ^ and ipj are independent and equally probable. We 
also consider that the sum (or difference) of two angles ip, as well as a product nip (where 
n is an arbitrary number) is also distributed on the interval (— 7r,7r). Another words, 
we deal with a wrapped angular distribution [20]. In this case, we can use two following 
arguments. First, the difference ip = ipi — ijjj of two uniformly distributed independent 
variables ipi and ipj also has a uniform distribution. In this case the probability density 
function (PDF) of ^ is a convolution of the corresponding PDFs of ipi and ipj [21], 
which is just a constant after wrapping to the interval ( — tt, 7r) . Second, using the rules 
of probability theory for transformation of random variables (2TJ |22], it can be shown 
that the product nip is also a uniformly distributed variable. Using a combination of 

= 0. 
N. 



these two arguments in (A.7) and (A. 8) we find that (C nii „ 2 ) M 
Therefore, the statistical estimate in (A.5) becomes a real number 



and (D„ 1>n2 ) 



■n-2/M 



Using this result in (A. 3) we obtain (13). 



(A nii „ 2 ) M 
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